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Description 



METHOD, SYSTEM AND STORAGE 

MEDIUM FOR REFERENCE 
NORMALIZATION FOR BLOCKED 
REFERENCE CHANNELS 

Background of the Invention 

[0001] This invention relates generally to a method and system 
for reference normalization in an x-ray image, and more 
particularly, to a method and system for reference nor- 
malization for blocked reference channels in an x-ray im- 
age generated by a computed tomography (CT) system. 

[0002] In CT imaging systems, an x-ray source projects a fan- 
shaped beam that is collimated to lie within an X-Y plane, 
generally referred to as an "imaging plane", of a Cartesian 
coordinate system toward an array of radiation detectors, 
wherein each radiation detector includes a detector ele- 
ment disposed within the CT system so as to receive this 
fan-shaped beam. An object, such as a patient, is dis- 



posed between the x-ray source and the radiation detec- 
tor array so as to lie within the imaging plane and so as to 
be subjected to the x-ray beam, which passes through the 
object. As the x-ray beam passes through the object, the 
x-ray beam becomes attenuated before impinging upon 
the array of radiation detectors. The intensity of the at- 
tenuated beam radiation received at the detector array is 
responsive to the attenuation of the x-ray beam by the 
object, wherein each detector element produces a sepa- 
rate electrical signal responsive to the beam intensity at 
the detector element location. These electrical signals are 
referred to as x-ray attenuation measurements or x-ray 
images. 

[0003] Moreover, the x-ray source and the detector array may be 
rotated, with a gantry within the imaging plane, around 
the object to be imaged so that the angle at which the x- 
ray beam intersects the object constantly changes. A 
group of x-ray attenuation measurements, i.e., projection 
data, from the detector array at one gantry angle is re- 
ferred to as a "view". A "scan" of the object comprises a 
set of views made at different gantry angles during one 
revolution of the x-ray source and the detector array. In 
an axial scan, the projection data is processed so as to 



construct an image that corresponds to a two-di- 
mensional slice taken through the object. In CT systems 
that employ a single detector array, the slice thickness is 
controlled and determined by the width of the collimator, 
while in CT systems that employ a multiple detector array, 
the slice thickness is controlled and determined by sum- 
ming the contributions of a plurality of detector sub-units 
and by physically moving the collimator to the outer 
edges of each slice. 
[0004] Reference normalization is an important step in CT pre- 
processing operations to reduce or remove the impact of 
the x-ray tube output fluctuation. For this purpose, a set 
of reference channels (also referred to herein as reference 
detectors) are placed slightly outside the reconstruction 
field of view (FOV) in the array of radiation detectors so 
that the reference channels receive x-ray photons directly 
from the x-ray source without interference from the 
scanned object. These reference channels monitor the x- 
ray tube flux and the measured signal is applied to the 
measured projections. The impact of any variations in the 
x-ray tube output on the measured projections is thereby 
properly removed. In a typical configuration, three refer- 
ence channels are placed at the left-hand side of the ar- 



ray. However, for larger than normal -si zed patients and/ 
or with the scanned object off-centered, there is a greater 
chance that these reference channels are blocked during 
the scan. In addition, with the introduction of wide-bore 
scanners that cover a larger FOV with the same detector, 
the reference channels are routinely blocl<ed by the pa- 
tients. To ensure the proper functioning of the reference 
normalization, different correction steps have to be devel- 
oped to account for situations where the reference chan- 
nels are blocked. 
Brief Description of the Invention 

[0005] One aspect of the invention is a method for performing 
reference normalization for a projection created by an 
imaging system. The imaging system includes a radiation 
detector array with a right and left edge. The method in- 
cludes receiving a projection dataset created by the imag- 
ing system in response to a varying x-ray tube current. 
The projection dataset includes a view and predicted 
fluxes are calculated for the set of reference channels 
within the view. A right set of reference channels are lo- 
cated proximate at the right edge of the detector array 
and a left set of reference channels are located proximate 
to the left edge of the detector array. The method also in- 



eludes calculating average actual fluxes for the sets of 
reference channels. Based on the predicted reference 
fluxes and the average fluxes, a reference correction value 
for the view is determined. The correction value is applied 
for the reference correction of the view. 
[0006] Another aspect of the invention is a method for perform- 
ing reference normalization for a projection created by an 
imaging system that includes a radiation detector array 
with a right and left edge. The method includes operating 
the imaging system so as to create a projection dataset 
responsive to an object. The imaging system is operated 
at varying x-ray tube currents to create the projection 
dataset. The projection dataset includes a view. Predicted 
fluxes for the sets of reference channels within the view 
are calculated, wherein a right set of reference channels is 
located proximate at the right edge of the detector array 
and a left set of reference channels is located proximate 
to the left edge of the detector array. Average actual 
fluxes are calculated for the sets of reference channels. A 
reference correction value is determined for the view 
based on the predicted reference fluxes and the average 
actual fluxes. The reference correction value is applied to 
the view. 



[0007] A further aspect of the invention is a medium encoded 

with a machine-readable computer program code for per- 
forming reference normalization for a projection created 
by an imaging system that includes a radiation detector 
array with a right and left edge. The medium includes in- 
structions to implement a method including receiving a 
projection dataset created by the imaging system in re- 
sponse to a varying x-ray tube current, where the projec- 
tion dataset includes a view. The predicted fluxes for the 
sets of reference channels within the view are calculated. 
A right set of reference channels is located proximate at 
the right edge of the detector array and a left set of refer- 
ence channels is located proximate to the left edge of the 
detector array. The average actual fluxes are calculated 
for the sets of reference channels. The referenced cor- 
rected values are then determined based on the predicted 
reference fluxes and the average actual fluxes. The refer- 
ence correction value is applied to the view. 

[0008] Another aspect of the invention is a system for performing 
reference normalization for an image. The system in- 
cludes a gantry having an x-ray source and a radiation 
detector array, wherein the gantry defines an object cav- 
ity. The x-ray source and the radiation detector array are 



rotatingly associated with the gantry so as to be separated 
by the object cavity. The detector array includes a right 
and left edge. An object support structure movingly asso- 
ciated with the gantry allows communication with the ob- 
ject cavity. The system also includes a processing device 
including instructions to implement a method that in- 
cludes receiving a projection dataset created by the imag- 
ing system in response to a varying x-ray tube current. 
The projection dataset includes a view. Predicted fluxes 
for the sets of reference channels are calculated within the 
view, where a right set of reference channels is located 
proximate at the right edge of the detector array and a 
left set of reference channels is located proximate at the 
left edge of the detector array. Average actual fluxes for 
the sets of reference channels are calculated. A reference 
correction value for the view is determined based on the 
predicted reference fluxes and the average actual fluxes. 
The reference correction value is applied to the view. 
[0009] Another aspect of the invention is a system for performing 
reference normalization for a projection. The system in- 
cludes an imaging system. The imaging system includes a 
radiation detector array with a right and left edge. The 
imaging system also includes an object disposed so as to 



be communicated with tlie imaging system. Tlie imaging 
system generates a projection dataset responsive to the 
object and to a varying x-ray tube current, where the pro- 
jection dataset includes a view. The processing device re- 
ceives the projection created by the imaging system. The 
processing device calculates predicted fluxes for the sets 
of reference channels within the view, where a right set of 
reference channels is located proximate at the right edge 
of the detector array and a left set of reference channels is 
located proximate at the left edge of the detector array. 
The processing device also calculates average actual 
fluxes for the sets of reference channels. A reference cor- 
rection value is determined for the view based on the pre- 
dicted reference fluxes and the average actual fluxes. The 
reference correction value is applied to the view. 
Brief Description of the Drawings 

[0010] Referring to the exemplary drawings wherein like ele- 
ments are numbered alike in the several FIGs.: 

[001 1] FIG. 1 is a perspective view of a CT imaging system and a 
patient disposed for imaging, in accordance with an ex- 
emplary embodiment of the present invention; 

[0012] FIG. 2 is a block schematic diagram of a CT imaging sys- 
tem, in accordance with an exemplary embodiment of the 



present invention; and 
[0013] FIG. 3 is a flow diagram of a metliod for reference nor- 
malization for blocked reference channels in a CT imaging 
system in accordance with an exemplary embodiment of 

the present invention. 
Detailed Description of the Invention 

[0014] An exemplary embodiment of the present invention in- 
cludes identifying a condition where the reference channel 
is blocked. When this condition occurs, the measured ref- 
erence channel, or reference detector, readings should not 
be used for the normalization process since they no 
longer represent the x-ray output flux. Some of the newer 
imaging systems, such as models of the wide-bore scan- 
ners include features where the x-ray tube current (and 
therefore the output flux) changes constantly so that the 
average projection noise can be held roughly constant 
when scanning objects that vary significantly. For exam- 
ple, when a scan is performed from the top of the head to 
the bottom of the feet of a patient, the x-ray tube current 
will vary as a function of where the current view is located 
within the anatomy (e.g., neck, stomach). In another ex- 
ample, the x-ray tube current will vary based on the angle 
of the patient relative to the x-ray source (e.g., side to 



side, front to back). Therefore, an independent signal tliat 
is closely related to the output x-ray flux and, at the same 
time, is not influenced by the patient blockage must be 
obtained. In an exemplary embodiment of the present in- 
vention, the sampled x-ray tube current is utilized as the 
independent signal. During CT scanning, the input current 
to the x-ray tube is sampled for every view and is stored 
in the scan file, or projection dataset. The projection 
dataset includes data corresponding to one or more views 
associated with a scan. The conversion factor that maps 
the tube current to the output flux is also stored in the 
file, since this factor is fixed for a known scanning proto- 
col (kV, slice thickness, etc.). Note that this factor varies 
between reference channels and CT systems because dif- 
ferent detectors have different gains to the x-ray input. 
[0015] Referring to FIG. 1 and FIG. 2 a representative CT imaging 
system 1 is shown and includes a gantry 2 having an x- 
ray source 4, a radiation detector array 6, a patient sup- 
port structure 8 and an object cavity 10. The x-ray source 
4 and the radiation detector array 6 are opposingly dis- 
posed so as to be separated by the object cavity 10. An 
object, such as a patient 12, is disposed upon the patient 
support structure 8 which is then disposed within the ob- 



ject cavity 10. The x-ray source 4 projects an x-ray beam 
14 toward the radiation detector array 6 so as to pass 
through the patient 12. The x-ray beam 14 is then colli- 
mated by a collimator so as to lie within an X-Y plane of a 
Cartesian coordinate system referred to as an "imaging 
plane". After passing through and becoming attenuated by 
the patient 12, the attenuated x-ray beam 16 is received 
by the radiation detector array 6. The radiation detector 
array 6 includes a plurality of detector elements 18. Each 
of the detector elements 18 receives the attenuated x-ray 
beam 16 and produces an electrical signal responsive to 
the intensity of the attenuated x-ray beam 16. In addition, 
some portions of the x-ray beam 16 may result in periph- 
ery sections of beams that don't pass through and be- 
come attenuated by the patient 12. These periphery sec- 
tions 46 may be received by reference channels, or refer- 
ence detectors 44, that are included in the radiation de- 
tector array 6 as depicted in FIG. 2 or by detector ele- 
ments 18 that are at the outer edge of the radiation de- 
tector array 6. 

[0016] In addition, the x-ray source 4 and the radiation detector 
array 6 are rotatingly disposed relative to the gantry 2 and 
to the patient support structure 8, so as to allow the x-ray 



source 4 and the radiation detector array 6 to rotate 
around the patient support structure 8 when the patient 
support structure 8 is disposed within the object cavity 
10. The x-ray projection data is obtained by rotating the 
x-ray source 4 and the radiation detector array 6 around 
the patient 12 during a scan. The x-ray source 4 and the 
radiation detector array 6 are in communication with a 
control mechanism 20 associated with the CT imaging 
system 1. Control mechanism 20 controls the rotation and 
operation of x-ray source 4 and radiation detector array 
6. 

[0017] jhe control mechanism 20 includes an x-ray controller 22 
that is in communication with the x-ray source 4, a gantry 
motor controller 24, and a data acquisition system (DAS) 
26 in communication with the radiation detector array 6. 
The x-ray controller 22 provides power and timing signals 
to the x-ray source 4 and the gantry motor controller 24 
controls the rotational speed and angular position of the 
x-ray source 4 and the radiation detector array 6. The 
DAS 26 receives the electrical signal data produced by the 
detector elements 18 and converts this data into digital 
signals for subsequent processing. The CT imaging sys- 
tem 1 also includes an image reconstruction device 28, a 



data storage device 30 and a processing device 32. Tlie 
processing device 32 is in communication witli tlie image 
reconstruction device 28, tlie gantry motor controller 24, 
the x-ray controller 22, the data storage device 30, an in- 
put device 34 and an output device 36. Moreover, the CT 
imaging system 1 also includes a table controller 38 in 
communication with the processing device 32 and the pa- 
tient support structure 8, so as to control the position of 
the patient support structure 8 relative to the object cavity 
10. 

[0018] The patient 12 is disposed on the patient support struc- 
ture 8, which is then positioned by an operator via the 
processing device 32 so as to be disposed within the ob- 
ject cavity 10. The gantry motor controller 24 is operated 
via the processing device 32 so as to cause the x-ray 
source 4 and the radiation detector array 6 to rotate rela- 
tive to the patient 12. The x-ray controller 22 is operated 
via the processing device 32 so as to cause the x-ray 
source 4 to emit and project a collimated x-ray beam 14 
toward the radiation detector array 6 and hence toward 
the patient 12. The x-ray beam 14 passes through the 
patient 12 so as to create an attenuated x-ray beam 16, 
which is received by the radiation detector array 6. 



[0019] The detector elements 18 receive the attenuated x-ray 
beam 16, produce electrical signal data responsive to the 
intensity of the attenuated x-ray beam 16 and communi- 
cate this electrical signal data to the DAS 26. The DAS 26 
then converts this electrical signal data to digital signals 
and communicates both the digital signals and the electri- 
cal signal data to the image reconstruction device 28, 
which performs high-speed image reconstruction. This 
information is then communicated to the processing de- 
vice 32, which stores the image in the data storage device 
30 and displays the digital signal as an image via the out- 
put device 36. 

[0020] FIG. 3 is a flow diagram of a method for reference nor- 
malization for blocked reference channels in a CT imaging 
system in accordance with an exemplary embodiment of 
the present invention. The process depicted in FIG. 3 is 
initiated for each view within a projection dataset. At step 
302, a predicted reference flux is calculated for the radia- 
tion detectors located at the leftmost side of the radiation 
detector array 6 (e.g., the reference channels), as well as 
the radiation detectors located at the rightmost side of the 
radiation detector array 6. The predicted x-ray flux ob- 
served at the reference channels, f (t), may be derived as 



follows: 

// (0 = 9- c(t) ; 

where c(t) is the measured x-ray tube current curve, q is 
the conversion factor and t is the view sampling index. 
Because the averaged reference readings are utilized for 
normalization, the detector channel index is dropped. The 
subscript "/ " in the above equation denotes the fact that 
most CT scanners include reference channels on the left 
hand side of the detector. 
[0021] The predicted flux, (t), is also calculated for the right- 
most channels. In an exemplary embodiment of the 
present invention, the three rightmost channels in the de- 
tector array are utilized as reference channels. It should 
be noted that (Dand (t) are related by a simple scaling 
that reflects different detector gains. This scaling factor 
may be determined by the ratio of the left and right chan- 
nel readings during air calibration. Experiments have 
shown that the relationship between the recorded tube in- 
put value and the measured reference value is not linear. 
That is, the value of q changes with the change of x-ray 
tube current. Therefore, q(c) is a tube input dependent 



function. The function can be determined experimentally 
or theoretically. For example, experiments have shown 
that the following function holds for a wide range of x-ray 
tube current: 

q(c) = qo[\-sic{t)-Co)]; 

where q the conversion factor q(c) at a particular, 
known x-ray tube current value (e.g., q ^ = 8.5), 5 is a 
scaling parameter (e.g., s = 1.0x10'^), and c ^ is a tube in- 
put reference value (e.g., c ^ = 9915). Alternatively s can 

be made a function of c .It should be noted that the con- 

0 

version factor q also changes with the DAS sampling rate. 
(This simple linear relationship will not be discussed 
here.) Note that here the tube input value is in some mea- 
surement unit and not the x-ray tube current milli-amp. 
The measurement unit can be converted to the x-ray tube 
current in mA by a simple scaling. 
[0022] Next, at step 304, the average measured flux, or average 
reference channel reading, for both the left side and the 
right side is calculated. The average left reference channel 



reading, R ^ (t), may be expressed as: 

m ^ 

i=\ 

where, p(i,t) is tlie detector reading for cell i and view t, 
and m is the parameter that controls the number of chan- 
nels to average. In an exemplary embodiment of the 
present invention, there are three reference channels 
(m=3). Note that in these equations, the detector row no- 
tation is dropped for simplicity. It is understood that this 
operation needs to be repeated for each row in the detec- 
tor. 

[0023] At this point, the predicted value could be compared 

against the measured value for reference blockage. How- 
ever, exemplary embodiments of the present invention 
perform additional steps to verify the data and to avoid 
cases where the reference channels are slightly blocked 
and non-optimal reference values are used for normaliza- 
tion. Although there is only one set of reference channels, 
or reference detectors 44, located on the left side of the 
detector, the rightmost detector channels are utilized as 



temporary reference channels. This is based on the obser- 
vation that unless the scanned object is very large, the left 
channel blockage is likely caused by the off-centered po- 
sitioning of the patient. Therefore, there is a good proba- 
bility that the right side channels are not blocked. 
[0024] The right reference reading, R (t), may be represented as: 



where n is a parameter that controls the number of chan- 
nels to average (e.g., n=3), and Nis the total number of 
detector channels in a row. Since the detector gains for 
the left reference channels and the right-most channels 
are likely to be different, R ^ (t) is scaled by a scaling fac- 
tor. This scaling, or gain, factor is determined by the ratio 
of the left and right channel readings during air scan (no 
object in the beam path). For simplicity, R ^ (t) is utilized to 
denote the properly scaled average reading of the right- 
most channel. 

[0025] At step 306, the predicted flux values (f (t)aLnd (Ware 




N 



n 



checked and corrected for errors in one or more of the 
conversion factors. In general, the predicted values based 
on the x-ray tube current, f ^(t) and ^ (t), are used as the 
bases of reference correction. However, to guard against 
the situation in which either the conversion factor, gain 
and/or the current value are wrong, it is necessary to 
make sure that the magnitude of these values is at least a 
significant fraction of the measured x-ray flux (R ^ (t) andR 
^ (t)) at the reference channels. To accomplish this, the fol- 
lowing operation is performed: 

f//(0, Mt)>aRi(t) m>aRr{t) 

In an exemplary embodiment of the present invention, the 
parameter a, is set to 0.95 and the scaling parameter e, is 
set to 0.97. It should be noted that f ^(l) and R ^ a; would 
be equal under non-blocked conditions. 
[0026] At step 308, the least blocked side is selected to be uti- 
lized as the reference correction value. The ratio of the 
measured reference signal against the predicted reference 
signal is calculated for both the left side and the right side 
to ensure that the side with the highest ratio is always 
used for reference correction. This is to ensure that if one 
of the two sides is blocked by the patient (even if it is only 



slightly), the unblocked side is used. The ratios, □ ^ (t)and 
□ (t), are calculated as follows: 



ratio is used for the reference correction. For illustration, 
assume that ri^(t)>r||(t). (The other case can be treated in 
exactly the same manner by substituting the subscript "r" 
by "r in all the following equations.) 
[0027] At step 310, the reference correction value for the previ- 
ous view scaled by the x-ray tube current change is uti- 
lized if both sides are blocked. This determination of the 
reference correction value, wft), can be expressed as the 
equation: 





ri|(t) is compared against r|^(t) and the side with a higher 




where 6 is a parameter (e.g., 6=0.95 in an exemplary em- 



bodiment of the present invention). Tiiis equation indi- 
cates that the measured reference value should be utilized 
only if it is more than 95% of the predicted signal based 
on the x-ray tube current. If this condition does not hold, 
it indicates that channels, or detectors, on both sides of 
the detector array 6 are blocked (since the higher of the 
two sides does not meet the requirement) and that the 
measured reference value is not reliable. In the case where 
both sides are blocked, the reference correction value of 
the previous view, w(t-l), scaled by the x-ray tube current 
change c(t)/c(t-l) is used as the reference correction value 
w(t). For the final reference correction, the measured pro- 
jection is divided by the correction value, w(t). 
[0028] In an exemplary embodiment of the present invention, the 
processing described above in reference to FIG. 3 may be 
implemented by computer instructions located on the 
processing device 32 with data being stored on the data 
storage device 30. In an alternate exemplary embodiment, 
the processing described in FIG. 3 may be implemented by 
a processor located remotely from the processing device 
32. 

[0029] The formulas described in reference to FIG. 3 are utilized 
by an exemplary embodiment of the present invention. As 



is known in tlie art, other formulas may be utilized to per- 
form the steps described in reference to FIG. 3. For exam- 
ple, in the w(t) calculation, only one preceding view w(t-l) 
is used when both sides are blocked. To improve robust- 
ness of the algorithm, multiple preceding views can be 
used. In the example listed above, only the left reference 
channel is designed into the detector. However, one can 
design a set of right reference channels to the right of the 
detector to improve the performance, since the further 
away a channel is from the center, the less is the proba- 
bility that the channel will be blocked. 

[0030] In accordance with an exemplary embodiment of the 

present invention, while a method and system for refer- 
ence normalization for blocked reference channels in an 
x-ray image, generated by an imaging system is de- 
scribed and discussed herein with reference to a com- 
puted tomography (CT) imaging system, it should be un- 
derstood that the method and system of the invention 
may be applied to other imaging systems (e.g.. Magnetic 
Resonance Imaging (MRI)). 

[0031] The ability to perform reference normalization for blocked 
reference channels when the x-ray tube current has the 
potential of changing from view to view can lead to more 



accurate images being produced by tlie imaging system. 
In addition, tlie ability to use regular channels, or detec- 
tors, that are located on the radiation detector array to 
supplement the reference channels may lead to less lost 
data due to reference channel blockage. Further, the abil- 
ity to utilize a previous reference correction value when 
channels at both sides of the radiation detector array are 
blocked may also lead to improved performance. 
[0032] As described above, the embodiments of the invention 
may be embodied in the form of computer-implemented 
processes and apparatuses for practicing those processes. 
Embodiments of the invention may also be embodied in 
the form of computer program code containing instruc- 
tions embodied in tangible media, such as floppy 
diskettes, CD-ROMs, hard drives, or any other computer- 
readable storage medium, wherein, when the computer 
program code is loaded into and executed by a computer, 
the computer becomes an apparatus for practicing the in- 
vention. An embodiment of the present invention can also 
be embodied in the form of computer program code, for 
example, whether stored in a storage medium, loaded into 
and/or executed by a computer, or transmitted over some 
transmission medium, such as over electrical wiring or ca- 



bling, through fiber optics, or via electromagnetic radia- 
tion, wherein, when the computer program code is loaded 
into and executed by a computer, the computer becomes 
an apparatus for practicing the invention. The technical 
effect of the executable instructions is to enhance the 
quality of the x-ray image when one or more of the refer- 
ence channels are blocked, resulting in less radiation ex- 
posure by decreasing the number of images taken. When 
implemented on a general-purpose microprocessor, the 
computer program code segments configure the micro- 
processor to create specific logic circuits. 
[0033] While the invention has been described with reference to 
exemplary embodiments, it will be understood by those 
skilled in the art that various changes may be made and 
equivalents may be substituted for elements thereof with- 
out departing from the scope of the invention. In addition, 
many modifications may be made to adapt a particular 
situation or material to the teachings of the invention 
without departing from the essential scope thereof. 
Therefore, it is intended that the invention not be limited 
to the particular embodiment disclosed as the best mode 
contemplated for carrying out this invention, but that the 
invention will include all embodiments falling within the 



scope of the appended claims. Moreover, the use of the 
terms first, second, etc. do not denote any order or im- 
portance, but rather the terms first, second, etc. are used 
to distinguish one element from another. 



